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Abstract. The blazar CTA 102 underwent a major radio flare in April 2006. We used 
several 15 GHz VLBI observations from the MOJAVE program to investigate the influence 
of this extreme event on jet kinematics. The result of modeling and analysis lead to the 
suggestion of an interaction between traveling and standing shocks 0.2 mas away from the 
VLBI core. 
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1. Introduction 

The radio source CTA 102 (z — 1 .037) is one of 
the most observed AGN in the northern sky. Its 
observational history sta rted in the late 1950s 
dHarris & Robertsl [i9 60). Flux density varia- 
tions were reported in the source dSholomitskiil 
119651) . which lead other authors to suggest that 
the signal w as coming from a n extraterrestrial 
civilization (lKardashevlll964 . Later, CTA 102 
was identified as a quasar. Due to its strong flux 
density variability CTA 102 has been the target 
for numerous observations at different wave- 
lengths. 

On kpc-scales CTA 102 consists of a centra l 
core and two faint lobes (Spen cer etal] [l989). 
The brighter lobe has a flux density of 170 mJy 
at a distance of 1 .6 arcsec from the core at po- 
sition angle (P. A.) of 143° (measured from 
North through East). The other lobe, with a flux 
density of 75 mJy, is located 1 arcsec from the 



center at P. A. -43°. The spectral indices of the 
lobes are -0.7 for the brighter one and -0.3 for 
the other. 

High resolution 15 GHz VLBI images show a 
curved jet with components which exhibit ap- 
parent velocities up to 15.4 + 0.9 c, adopting a 
cosmology with O,,, = 0. 27, Qa = - 73 an d 
H = 71kms~ 1 M pc~ 1 dLister et all 12009 ' 



iJorstad et all d2005l) and iHovatta et all d200! 
derived bulk Lorentz factors, F, between 15 
and 17 and Doppler factors, 5, between 15 and 
22 from 43 GHz VLBI observations and from 
single-dish light curves. 
A major flux density outburst occurred in 
April 2006, which offers a unique oppor- 
tunity t o study CTA 102 u nder these con- 
ditions ( From met al. I l2010l) . For our anal- 
ysis we used data from the monitoring of 
CTA 102, performed within the framework of 
the 15 GHz MOJAVrfl program (Monitoring 
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of Jets in A ctive galactic nucle i with VLBA 
Experiments: llister et al.1 d2009l) '). 

2. Observations 

The MOJAVE observations of CTA 102 be- 
tween 2005 and 2008 have been used to an- 
alyze kinematic changes in the source during 
the 2006 radio flare. The raw data were cal- 
ibrated using standard AIPS procedures. The 
calibrated data were fitted by several circular 
Gaussian components using DIFMAP. The fit- 
ted components are labeled from CI to CI 2, 
in inverse order of distance to the core. In 
Fig. [1] one can see a high resolution image of 
CTA 102 observed at 15 GHz on 6th of January 
2007. 
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Fig. 1.15 GHz VLBI image of CTA 102 with fitted 
circular Gaussian components observed on 6th of 
January 2007. The map peak flux is 1.13 Jy /beam, 
where the convolving beam is 1.1 x 0.6 mas at P. A. 
-4.8°. 

Using the evolution of the flux density and 
displacement of the fitted components from the 
VLBI core we could identify the modeled fea- 
tures between the different epochs. This ap- 
proach allows us to identify the feature labeled 
C12 as a possible candidate for a newly ejected 



component. The evolution of the flux density 
and position of component C12 are presented 
in Fig. [2] 
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Fig. 2. Top panel: Evolution of the separation from 
the VLBI core of component C12. Bottom panel: 
Flux density evolution of component C12, the VLBI 
core and the total flux density. The uncertainties in 
the flux densities are around 5% and the positional 
errors are ~ 1/5 of the restoring beam size (typically 
around 0.01 mas) 

The separation from the core of component 
C12 remains constant until 2006.6, followed 
by a sharp acceleration an d strong decrease in 
the observed flux density. iLister et all d2009t) 
derived an apparent speed of 8.5 ± 1.0 c and 
an ejection time of f e j = 2005.04 + 0.3 yr in- 
cluding the stationary phase before mid 2006 
in the fit. For our analysis we separated the C12 
trajectory into a stationary section, C12 a , (be- 
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fore 2006.6) and moving section, C12b (after 
2006.6). For CI 2b we derived an apparent ve- 
locity of 17.3 ± 0.7 c and an ejection time of 
<ej =2006.16 ±0.05 yr. 

3. Discussion 

There is a flare in 15 GHz observations of 
CTA 102 around 2005.8 (see increase in to- 
tal and core flux density in bottom panel of 
Fig.|2]i. Furthermore, as already mentioned, the 
trajectory of component C12 is stationary until 
2006.6. In the evolution of the flux density of 
this component one can see a global maximum 
around this time (see Fig . |2j . 
A possible scenario which explains the trajec- 
tory of C12 and the evolution of the flux den- 
sity could be an interaction between traveling 
and standing shocks. 

The observed flare in CTA 102 could be ex- 
plained by the propagation of a relativistic 
shock wave. This wave is generated by pres- 
sure mismatches at the jet nozzle and during 
its way downstream it reaccelerates the under- 
lying jet particles at the shock front. This inter- 
action of the shock with the jet flow leads to an 
increased emissivity and can be ob served with 
VLBI techniques (lMarscheJ2009h . 
The stationarity of component C12 a could 
be associated with a standing shock around 
0.15 mas away from the core (see Fig. [2]). A 
standing shock is a stationary wave in non- 
pressure matche d jets, created during the colli- 
mation of the jet dDalv & Marscherll988l:lFallel 
1991). The position of such a shock remains 
stationary as long as the relative conditions 
of the jet and ambient medium remain un- 
changed. 

Using the trajectory of component CI 2b, we 
derived an ejection time of f e j = 2006.16 + 
0.05 yr. From this result it is not immediately 
obvious that one can associate the CI 2b with a 
travelling shock wave which triggered the ob- 
served flare around 2005.8. This temporal dis- 
crepancy between the onset of the flare and the 
calculated ejection time of component C12b 
could be explained in the following way: Due 
to the resolution limit of the 15 GHz VLBI 
observations there is a lack of data less than 
0.15 mas from the core. By using a linear in- 



terpolation based on the observed trajectory of 
component CI 2b, we implicitly assumed a con- 
stant velocity for this compone nt and do not 
take ac celeration into account dHoman et al. I 
20090). If there is acceleration in the trajectory 
of C12b, this feature could have been ejected 
earlier than the calculated 2006.16. Therefore 
we could identify the CI 2b trajectory with the 
interaction of a traveling shock wave with the 
underlying jet flow causing the 2005.8 flare. 
Such a moving shock wave would reach its 
peak flux density soon after ejection from the 
core, followed by a permanent flux density de- 
crease due to, e.g., adiabatic expansion. In the 
case of CTA 102 there is a stationary feature 
around 0.15 mas from the core and therefore 
this effect can not be observed due to blend- 
ing with the emission coming form the station- 
ary feature (nearly constant flux density after 
2006.2 of component C12). While the travel- 
ing shock is propagating downstream, it will 
collide with the stationary one after some time 
t co i = 2006.6. During this collision the trav- 
eling shock will encounter a local increase in 
the particle density compared to the surround- 
ing flow. The higher density at this position 
would then lead to an increase in the emis- 
sivity because of the shock front accelerating 
more particles. This process could explain the 
global flux density maximum of component 
C12 around 2006.6. Furthermore, the stand- 
ing shock would be pulled away by the mov- 
ing one and after some time the stationary 
feature should re-ap pear at the same position 
dMimica et al.ll2009h . 

This scenario is in good agreement with the ob- 
served flux density evolution of the component 
C12. We could not detect the stationary sec- 
tion, C12 a , after the collision of the two waves 
in 2006.6. A possible explanation could be the 
limited dynamic range of the VLBI observa- 
tions or changes in the ambient medium con- 
figurations during the collision of the traveling 
shock wave. 

To confirm our suggestion of a shock-shock in- 
teraction in the jet of CTA 102 one should per- 
form a spectral analysis on the available single- 
dish and multi-frequency VLBI observations. 
The overall behaviour in the turnover fre- 
quency turnover flux density plane is in agree- 
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ment with the standard shock-in-jet model with 
a double hump feature. This feature could be 
explained by an increase in the Doppler fac- 
tor or by a shock-shock interaction. The results 
of this analysis will be published elsewhere. 
Future research on the shock-shock interaction 
should also include relativistic (magneto) hy- 
drodynamic (RMHD) simulations. These sim- 
ulations could help us to better understand 
the process during the interaction of the trav- 
eling shock wave with the stationary one 
and could bridge the observational time gaps. 
Furthermore these calculations could provide 
an estimate for the time of re-appearance of the 
standing shock wave. 

4. Conclusions 

We presented the scenario of an interaction be- 
tween traveling and standing shock as an ex- 
planation for the observed 2006 radio flare in 
CTA 102. Furthermore this scenario could ex- 
plain the trajectory of component C12, which 
could be associated with the flaring event. 
The combination of multi-frequency observa- 
tions (single-dish as well as multi-frequency 
VLBI observations) and R(M)HD simula- 
tions will allow us to confirm, or not, the 
assumption of a shock-shock interaction in 
CTA 102. The analysis of single-dish obser- 
vations in_Jhe_Jrameworkof_the shock-in-jet 
model dMars cher & Gearlll985l) will also pro- 
vide estimates of the evolution of the magnetic 
field, Doppler factor and spectral index of the 
relativistic electron distribution (Fromm et al. 
2010, in prep.). 
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